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where f{x) is the decreasing function that satisfies — = — x^, 
for < a < 6. When a is an integer and 6 = a + 1 we deduce several 
combinatorial results. These include an asymptotic formula for the 
number of integer partitions not having a consecutive parts, and a 
formula for the metastability thresholds of a class of threshold growth 
cellular automaton models related to bootstrap percolation. 

1 Introduction 

Let < a < b and define / = fa,b ■ [0, 1] — * [0, 1] to be the decreasing 
function that satisfies 

[/(x)r - = x» - x^ o<x<i. (1) 
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Figure 1: An illustration of the definition of the function /. 



y = x'^ — 




P f{u) 1 



X 



Our central result is the following. 
Theorem 1 For every < a < b, 



^ - log/(x) , TT- 



-dx 



X 



Sab 



Note that — x'' is increasing on [0, p] and decreasing on [p, 1] where 
P = Pa,b G (0, 1) is defined by 



-logp 



log b — log a 
b — a 



(2) 



It follows that / is uniquely determined by the above conditions, and satisfies 
/(O) = 1, /(I) = 0, f{p) = p, 

and 

f{f{x)) = x, 0<x<l. 

See Figure m 

When a is a positive integer and b = a + 1, Theorem ^ has the following 
consequences. 
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Probabilistic application. Let < s < 1, and let Ci, C2, ... be indepen- 
dent events with probabilities 



under a probability measure Pg. (We can think of Cn as the event that at 
least one occurs of a further set of n independent events each of probability 
s). Let /c be a positive integer, and let be the event 

Ak = f]{Ci u Ci+i u ■ ■ ■ u a+k-i) 
1=1 

that there is no sequence of k consecutive Cj's that do not occur. 
Theorem 2 For every positive integer k, 

The next two applications are consequences of Theorem |21 

Number-theoretic application. A partition of a positive integer n is 
an unordered multiset of positive integers (called parts) whose sum is n. Let 
Pkin) be the number of partitions of n that do not include any set of k distinct 
consecutive parts. (So for example P2(4) = 4, since the relevant partitions of 
4 are (4), (3, 1), (2, 2), (1, 1, 1, 1); the partition (2, 1, 1) is not allowed because 
it has the 2 consecutive parts 1,2). 

Theorem 3 For every integer k > 2, 



logpfc(n)~7r^/^(^l-^^^^^^^n as 00. 

Application to cellular automata. Threshold growth models are a class 
of simple cellular automaton models for nucleation and growth; see |2],|H1, 
in]; ^2] and the references therein. Elements of the two-dimensional integer 
lattice are called sites. At each time step t = 0, 1, 2, . . ., a site is either 
active or inactive. A site z has a neighborhood N{z) C defined by 

N{z) = {z + w: w e N}, 
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where A^(= A^(0)) is some fixed finite subset of Tl? . We also fix an integer 
Q called the threshold. The system evolves over time according to the 
following rules. 

(i) A site that is active at time t remains active at time t + 1. 

(ii) A site z that is inactive at time t becomes active at time t + 1 if and 
only if its neighborhood iV(2;) contains at least Q active sites at time t. 

Consider the random initial state in which at time 0, each site in the 
L hj L square is active with probability s, independently for 

different sites, while all sites outside the square are inactive. Let I{L, s) 
be the probability that every site in the square eventually becomes active. 
A central question is to determine for various models the behavior of the 
function I{L,s) as L oo and s — simultaneously. 

Theorem 4 Let k > 2 be an integer, and consider the threshold growth 
model with neighborhoods given by 



and threshold 6 = k. For L ^ oo and s ^ simultaneously we have 

(i) z/liminf slogL > A then I{L, s) ^ 1; 

(ii) if limsup s logL < A then I{L, s) 0, 



Further integrals. We can evaluate several other definite integrals using 
Theorem m Recall the definition of p in Q. 

Theorem 5 For every < a < b, 



N = Nk = {{-V, 0), {v, 0), (0, -v), (0, v):v = l,2,... 



k-1} 



where 



A 



3k{k + l)' 




and 
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Define / : [0, 1] [0, 1] to be the decreasing function tliat satisfies 

— /(x)log/(x) = —xlogx, < X < 1. (3) 



Theorem 6 



Theorem 7 







^ - log f{x) 

X 6 



log/(x)^_^ TT^ 1 



and 



X 6 2' 



6 2 



Remarks. As we shall see in the proof of Theorem d the result for fa^b 
implies that for fa-y,!)^ for any 7 > via an easy argument. The case a = 
1,6 = 2 is easy; in that case we have f{x) = 1 — x, and the integral in 
Theorem^Jis standard (j7j, number 4.291.2). The case a = 2,6 = 3 also has 
an explicit formula for /, and this was used to prove Theorem C] in that case 
in ^21- Our proof of the general case uses an entirely different approach. 

The case /c = 2 of Theorem IHl can also be deduced from a certain partition 
identity (see Section |^. This raises the possibility of a family of partition 
identities corresponding to other values of k. If such identities could be found, 
they might also lead to alternative (combinatorial) proofs of Theorems I1I2I 
We discuss these matters in more detail in Section EJ 

The case /c = 2 of the threshold growth model in Theorem 0] is called 
bootstrap percolation. Theorem |3] was proved for that case in The 
general version is proved by a modification of the proof in ^21; making use 
of Theorem 12 above to obtain the numerical value of A. In Section El we give 
an account of the proof, omitting some of the details. 

Prior to the proof in even the existence of the sharp constant A in 
Theorem |3] was not known. On the other hand, analogues of Theorem |3] 
with two different constants Ai, A2 in (i),(ii) were known for a wide class of 
models. In some cases, the "scaling function" s log L is replaced by a different 
function of s, L. In particular, two-dimensional models are studied in detail 
in [HI, El- For neighborhoods as in Theorem 0] and threshold k < 6 < 2k — 2 
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for example, the results in those articles imply that the appropriate scaling 
function is s^"'^"'"^ logL. (The cases 9 < k and 9 > 2k — 2 turn out to 
be less interesting; in the former case, an active square of side k will grow 
forever, while in the latter case an inactive square of side k will remain 
inactive forever). The reason for the particular choices of A^, 9 in Theorem |3] 
is that our methods (combined with those of ^2]) yield the sharp constant A 
relatively easily in these cases. The extension to other A^, 9 remains an open 
problem. 



2 Integrals 



In this section we prove Theorems | 

It suffices to prove the Theorem ^ for the case b = a + 1. To check this, 
suppose that it holds for a given choice of a, b. For 7 > 0, let 

Replacing x with x"^ in ((T]), we see that 

g^^ix) - g'^ix) = x"^ - x^, 

and g is decreasing, so g = fa-t,by Supposing that the theorem is true for 
fa,hi we will check it for g: 



X 7 Jo X 



-dy 



l"^ Jo y 'j'^Sab 3 (07) (67)' 

In the second step above, we have made the change of variable y = x"' . So 
we may without loss of generality take 

b = a + l. 

Note that in this case, p = a/b. We will use r(-) to denote the usual gamma 
function. 

The proof of Theorem Q is based on properties of the function 

Fix)-V^^^^^^ (4) 
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By Stirling's formula, 



T{be) 

r{a£)£\ 




2tiM 



a 



£] oo. 



Since the maximum value of x"' — = x"'{\ — x) on [0, 1] is 



a' 



a 



6"' 



the series in (0]) converges uniformly on [0, 1], and hence defines a continuous 
function there. Note that the same cannot be said for the series for F' . In 
fact, F' is not continuous at x = p; one can show using the proposition 
below that F'{p—) = 1/p and F'{p+) = —l/p. This singularity will play 
an important role in the analysis. The following result contains the main 
properties of F that will be needed in the proof of Theorem ^ 

Proposition 8 Let b = a + 1. The function F has the following properties. 

(i) F{f{x))=F{x) on [0,1]. 
(a) F{x) = — log/(x) on [0,p]. 
(Hi) F{x) = — logx on [p, 1]. 



Proof. Part (i) is immediate from and the fact that the series in 
depends on x through the expression — x''. Part (ii) is a consequence 

of (i) and (iii). To see this, take x G [0,p]. Then f{x) G [p, !]• By (iii), 

F{f{x)) = — log/(x). Now use (i). 
Turning to the proof of (iii), define 



(iv) 






for complex z in the connected component Q of the set 
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that contains the segment (p, 1]. Note that Q is contained in the right half 
plane since for 'Re{z) = p, — z\ > p'^(l — p) = a°'/b''. Therefore z"" can 

be defined unambiguously on Q as an analytic function that takes the value 
1 at z = 1, and F is then analytic in Q. The function — logz is also analytic 
in Q, and can be chosen to take the value at z = 1. So, it suffices to show 
that F{z) = —logz in a complex neighborhood oi z = 1. Write w = 1 — z, 
and consider the neighborhood of 



In N, 

-log(l - w) = V — . 



m 

m=l 



Also, in A^, the following rearrangement is justified by absolute convergence 
of the series involved: 



oo 



F{l-w) = 



e=i 



TjM) ((l-w)"w 
T(ai)£\ a£ 



j^T{ai)i\aif^^T{ai-k + l)k^^ ' 



where 



T{bi) 1 ^ Tiai+1) ^ 

r(M)(-i)'"-' 



oo 
m=l 



-'m 



:^£!(m-£)!r(fe^-m + l) 

So, it suffices to prove that 6^ = for m > 1. 

To do so, use the property r(a + 1) = aV{a) to rewrite as 



1=1 



The summand above that would correspond to £ = is — 1/m. Therefore, 
6m = l/"*^ is equivalent to 



m— 1 



£=0 ^ ' i=\ 
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Now write 



(1 - x^) 



b\m 



X 



E 



m 



-Ifx 



To check dU)), it is then enough to show that 



dx 



m—l 



X 



x=l 



Let 



1 / l 

h{x) = -(' 



x\ 1 — X 



X ^ 1 



and h{l) = fe™, so that 



(1 - x '') 

X 



b\m 



h{x){i-xy 



Since h is C°° in a neighborhood of x = 1, 



dx''' 



—h{x){l — x) 







x=l 



as required. 

To prove part (iv) of the proposition, we use the standard beta integral 

^ ai-i/. xa2-i , r(ai)r(a2) 
u [l — u) du 



(see [TT] pl48). This gives 
F{x) 



-dx 



X 



T{ai + 02) 



r(6£) [x''{l-x)Xl 
ai X 



S r(a£)£! 

^ r{bi)r{ai)r{£ + i) 



□ 
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Proof of Theorem [H As remarked at the beginning of the section, we 
may assume b = a + 1. By Proposition |H1 (iv), it suffices to show that 

'- log /M,,,, (7) 

X In X 



To do so, let 







X „ X 



X L X 



p 



By Proposition IHl (ii), 

h = Ji. (8) 
Making the substitution y = f{x) and then integrating by parts gives 

h = 1^ ^""^yj^dy = -(logp)' + ^2, (9) 

since the boundary term at y = 1 vanishes as a result of 

/(y)~(l-y)^/^ yU- 
By Proposition |H1 (iii) and another integration by parts, 

J, = =^dx = (10) 
Jp X 2 

Combining Q and (fTUI) gives h — h = 2J2, and this, together with gives 
+ /2 = 2Ji + 2J2, which is ©. □ 

Proof of Theorem \E[ This follows immediately from above (which 
holds for all < a < 6) and Theorem ^ D 

Proof of Theorems EHZl We will take the limit a/6 ^ 1 in Theorems 
11151 Let e G (0, 1), and write /(x) = /i_i+e(a^) and 



X — x^"^^ 
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so that by we have 

^{fix)) = ^(x). (11) 
Recall that (f{x) is increasing on [0,p] and decreasing on [p, 1]. Note that 

ip{x) —a; log a; as e | 

uniformly in x G [0, 1], and hence 

fix) fix) as e i 0. 

Note also (from (0)) that p | as e J, 0. 

We will use dominated convergence. First observe (by differentiating) 
that fix) is decreasing in e for each x. Hence 

X — x'^ < (fix) < —X log X. 

Therefore there exists a fixed constant c satisfying 

< c < < p < 1/2 < 1 - c < 1 
such that for all e G (0, 1) we have 

u/2 < ifiu) < ^/u for u G [0, c], 

and 

(1 - n)/2 < (^(n) < 1 - n for n G [1 - c, 1]. 

It follows from ()11|) and the definition of / that there exist fixed positive 
constants c', c" with c' < 1/4 such that for all e G (0, 1) we have 

{1 - 2y/x, X < c' 

c", c' <x<l-c' 

{'-f^)\ x>l-c'. 

Here, the bounds in the first and third cases are obtained by using the bounds 
for if above, and solving (1 — /(x))/2 < y/x and (1 — x)/2 < a/ fix) respec- 
tively. Therefore for all e G (0, 1) we have 

-log(l-2v/x) X < d 



x 



X > 1 — c', 



-21og((l-x)/2) 
11 



where d" > 0. The function on the right is integrable on [0, 1] (see [7] for 
example). Hence taking e | and using the dominated convergence theorem, 
Theorem IHl follows from Theorem Q and Theorem [7| follows from Theorem 
El □ 

One may ask what Theorem Q yields in the limit when a/6 — * (or oo). 
In fact, it yields nothing new. Taking a — >■ with ab = 1 say, an argument 
similar to the above shows that the limit of the integral is 7r^/3, using only 
the (easy) case a = 1, 6 = 2 of Theorem [T] 



3 Probability 

In this section we prove Theorem |21 

We say that a (finite or infinite) sequence of events has a fc-gap if there 
are k consecutive events none of which occur. Thus is the event that the 
sequence Ci, C2, . . . has no fc-gaps. 

Lemma 9 Let Wi, . . . ,Wn be independent events each of probability u e 
(0,1). Then the probability gn{u) that the sequence Wi, . . . ,Wn has no k- 
gaps satisfies 

[fk,k+i{i - u)r < 9n{u) < [/fc,,+i(i - u)r-'+\ 

Proof. Writing / = fk,k+i{x), we have from ((T)) that 
and rearranging gives 

Provided x 7^ p we have / 7^ x, so we may divide through by / — x to obtain 

= {l- x){f-^ + xf~^ + x^f^'' + ■■■ + x^-^), (12) 

and ()12p holds when x = p also by continuity (or Q). 

We now prove the statement of the lemma by induction on n. For n = 
0, . . . , — 1 we have (?„(«) = 1, so the statement holds because /fc,fc+i(l — m) G 
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(0, 1). For n > 0, we may compute gn+k by conditioning on the first Wi to 
occur: 

gn+k = ugn+k^i + (1 - u)ugn+k^2 + (1 - ufugn+k-'i H H (1 - uf~^ugn 

= (1 - x){gn+k-i + xgn+k-2 + x'^gn+k-3 H ^- x'^'^gn), 

where we have written x = 1 — u. Comparing this with we deduce that 
if the lemma holds for gn, ■ ■ ■ , gn+k-i then it holds for gn+k- D 

Proof of Theorem [21 The idea of the proof is that when s is small, 
Ps(C*n) varies only slowly with ra, so we may use Lemma IHl to deduce that 
Ps(y4fc) behaves approximately like H^i /fc,fc+i(l — PslC'n)), and this in turn 
may be approximated using the integral in Theorem Q (after a change of 
variable) . 

It is convenient to write 

q = -log(l - s) 

so that 'PsiCn) = 1 — e""'' and g ~ s as s — ^ 0. Note that the indicator of Ak 
is an increasing function of the indicators of Ci, C2, . . ., so if we increase (re- 
spectively decrease) the probabilities of the Ci while retaining independence 
then we increase (respectively decrease) the probability of A^. We write 

and let C^, C~ be independent events with probabilities 

P,(C+) = l-e-*"« {i-l)r<n<ir- 

, . ^ \ s < n < r, 

^s\^n) ^i_e-»r-<7 ir<n<{i+l)r, 

for z = 1, 2, . . .. Then we have Vs{C~) < P,(C„) < Ps(C+), and so Ps{A~) < 
Ps(^fc) < ); where (respectively A^) is the event that the sequence 

C^, C2, ■ ■ ■ (respectively Cf , Cg", . . .) has no A;-gaps. 

Now we may bound Ps{A^) above by the probability that 

for every i>l, C'(t-i)r+i5 ■ ■ - ^^i^ has no fc-gaps. 

And we may bound Ps(v4^) below by the probability that 

, . . . , C~ all occur, and 
for every i > 1, C~_^^, . . . , C~_^_^_;^ has no /c-gaps, and C:^_^_^ occurs. 
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Hence, applying Lemma El and writing / = fk,k+i we have 

oo oo 
i=l i=l 

hence 

oo 

_ A; + 1) ^ - log f{e-'^^) < - log PsiAk) 

i=l 

oo oo 

<_rlogs + J2- log(l - e''"^') + (r - 1) 5^ - log f{e-'^'^). (13) 



i=l 



Applying the change of variable x = e ^ to the integral in Theorem ^ 
(with a = k, b = k + 1) gives 



log/(e ')dz 



3k{k + l) 
and in the special case k = 1, 



foo 2 

-log{l- e-ndz = -. 



Using the fact that — log /(e ^) and — log(l — e ^) are decreasing functions 
of z, (fT^ implies 



r — k + 1 
rq 



POO 

/ -log/(e-^)d^<-logP,(A,) (14) 

J rn 



' rq 

I roc ^ J. _i foo 

<-rlogsH / -log(l-e ^)dz ^ / -log/(e ""^dz. 

rq Jo rq Jq 

Now we let s 0. Using the facts that g ~ s, r ~ s"-*^/^, and both integrals 
are convergent, we obtain that the upper and lower bounds in (|14j) are both 
asymptotic to 

I _ TT^ 1 

-log/(e ')dz 



?>k{k + l)s' 

and hence the same holds for — logR5(Afc). □ 
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4 Partitions 



In this section we prove Theorem El 

Lemma 10 For any k>2, Pk{n) is a non- decreasing function of n. 

Proof. For /c > 3, the following transformation defines an injection of 
the set of partitions of n not containing k consecutive parts into the set of 
partitions of n + 1 not containing k consecutive parts, thus establishing the 
claim. If the partition does not contain all of the numbers 2, 3, . . . , /c as parts, 
then we may add another part equal to 1, transforming the partition of n into 
a partition of n+ 1. If the partition does contain 2, 3, . . . , /c as parts, then we 
may transform it into a partition of n + 1 by taking one of the parts equal 
to 2 and changing it into a 3. It is easy to verify that this is an injection. 

It remains to prove the claim when k = 2. For that case, we define 
the following transformation taking partitions of n without two consecutive 
parts injectively into partitions of n + 1 without two consecutive parts. If 
the partition does not contain any 2's, then we may add a 1. If the partition 
does contain 2's, we add 3 to the largest part in the partition and remove 
one 2. (This fails for the special partition 2 = 2 of ri = 2, for that case verify 
the claim directly). □ 

Proof of Theorem El Denote by 



the generating function oi pk{n) {k fixed). Let p{n) be the total number of 
(unrestricted) partitions of n, and denote its generating function 



oo 



X' 



n=0 



oo 



G{x) = 5]]p(n)a;". 



n=0 



By IE], pl8, we have 




oo 



< X < L 



(15) 
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We now observe that Gk{x) is closely related to the probability Ps{Ak) in 
Theorem |21 Let s = 1 — x. We may write the event as a disjoint union 
over the countable set S of all binary strings 01020304 ■ ■ ■ G {0,1}^ that 
contain only finitely many O's, and in which there are never k consecutive 
O's, of the event 

i : ai=l i : ai=0 

(By the Borel-Cantelli lemma, with probability one only finitely many of the 
Cj's will fail to occur). Therefore 



0.10.213 



E p.f n n c{] 

203---65 \t : ai = l I : 0^=0 / 

n d-^') n 



E 

oia2a3---£5 



— X ) X 
i : ai=l i : 0^=0 



i=l 
1 



E n T 

aia2a3---ScS i : 0^=0 



X 



G{x) 

Gu{x) 
G{x) 



j2 n ix'+x''+x''+ 



(by (HID) 



aia2a3---65 i : 0^=0 



since on expanding out the sum of the products, the different ways to get 
correspond exactly to partitions of n without k consecutive parts (choosing 
the power of corresponds to choosing the number of times the part i 
appears in the partition). Now using Theorem |21 and the standard fact 
pl9) 



logG'(x) 



TT 



6{l-x)' 



as X t 1, 



we obtain 



logG'fc(x) 



k{k + l) J 1- X 



as a; T 1. 



We now use (a special case of) the Hardy- Ramanuj an Tauberian Theorem 
[TUj . which says that if H{x) = Yl'^=o ^nX"', where 6„ a positive non- decreasing 
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sequence, and log if (x) ~ c/(l — x) as x t 1, then log6„ ~ as n — oo. 

Theorem 01 follows (using Lemma El)- D 



The case k = 2 and partition identities. The special case k = 2 oi 
Theorem IHl can be deduced (and from it the corresponding cases of Theo- 
rems |21 and using the following elementary partition identity due to P. A. 
MacMahon (0 pl4, examples 9, 10). 

The number of partitions of n not containing I's and not contain- 
ing two consecutive parts is equal to the number of partitions of 
n into parts all of which are divisible by 2 or 3. 

Denote by r(n) the number of such partitions of n. It is straightforward 
to check that r(n) < r{n + 2) for all n, since given a no-ones, no-consecutive- 
parts partition of n one may add 2 to its largest part to turn it (injectively) 
into such a partition of n-|-2. Furthermore, we will argue that the restriction 
of containing no I's does not influence the exponential rate of growth of the 
partition counting function, since we have the inequalities 



For the non-obvious part r{n — 1) < P2{'n) of the lower bound, use the 
(injective) transformation that adds 1 to the partition if there are no 2's, and 
otherwise takes a 2 and adds it, together with an additional 1, to the largest 
part. For the upper bound, use the transformation that takes a partition and 
deletes all the I's. 

Let R{x) = be the generating function of r(n). By the above 

partition identity we have for < x < 1 



n 



max 



|r(n — l),r(n)} < P2{n) < ^^r(£). 



(16) 



oo 



1 



R{x) = n 



(1 -X6^+2)(1 -x6'^+3)(l -x6^+4)(l -x6*+6)' 



i=0 



or 



oo 



logi?(x)= 




i=2,3,4,6 i=0 
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It can be shown in a manner analogous to the asymptotic behavior of G{x) 
cited above that for any j > 1, 



oo 



1 




6 ' 6(1 



x) 



as X t 1. 



i=0 



Therefore, summing over j = 2, 3, 4, 6 gives 



log R{x) 



4 



as X t 1- 



6 ' 6(1 -x) 



By the Hardy-Ramanujan Tauberian Theorem apphed to the function R{x) + 
xR{x) (the generating function of the non- decreasing sequence r(n — 1) + 
r(n)), we get 



which by (fTBj) gives Theorem El for k = 2. 

Now Theorem 121 may be deduced by following the arguments of the proof 
of Theorem IHl in the opposite direction (with the slight adjustment of re- 
placing Ps{A2) with Ps(Ci n A2) to account for the modified definition of 
the partitions), and Theorem ^ may be deduced by following the steps of 
the proof of Theorem |21 (adapted to fit the modified statement) in the op- 
posite direction. Note also that, as a consequence of MacMahon's identity, 
Ps(Ci n A2) has the intriguing factorization 



P,(CinA2) = JJ(l-x^'^+^)(l-x^'=+^) = (l-x)(l-x^)(l-x^)(l-x^^)--- 



(where again x = 1 — s). Can this fact be given a direct probabilistic proof? 

We remark finally that, in light of the above argument and the neat form 
of the exponential growth constant for Pk{n) in Theorem |2l it is tempting to 
conjecture the existence of partition identities for other integer values of k 
that would give an alternative proof of Theorem El (and therefore of Theorems 
121 and for integer a = k and b = a + 1. This would imply, by analytic 
continuation, the general case of Theorem Q thus giving an independent 
proof of Theorem ^ Presumably, such partition identities would equate the 
number of partitions of n not containing k consecutive parts and possibly 



log (r(n — 1) + r{n 



)) ~ '^-\fn as n — * 00 
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satisfying some other "mild" conditions, with the number of partitions of 
n whose parts satisfy some congruence restrictions modulo k{k + 1) (there 
should be two forbidden congruence classes) and other mild restrictions. The 
discovery of such identities would be an interesting positive use of partition 
asymptotics in the study of partition identities. See [H] for an example of a 
negative use of partition asymptotics, where they were used to prove the non- 
existence of certain partition identities. Also, see |3], jS], [ISl for discussion 
of connections between partition theory and various models in geometric 
probability, random matrix theory and statistical mechanics. 

5 Cellular Automata 

In this section we describe the proof of Theorem HJ The argument is a modi- 
fication of that in ^2] • We therefore omit many of the details, concentrating 
instead on the differences compared with ^21- The basic strategy is as fol- 
lows. Roughly, an L by L square will become fully active if and only if 
it contains at least one "nucleation center". A nucleation center should be 
thought of as a local configuration of active sites that will grow to occupy 
the entire square. One way for a nucleation center to occur involves the oc- 
currence of (two independent copies of) the event Ak in Theorem |21 Hence 
the estimate of Ps{Ak) in Theorem [pleads to the value of A. The above ideas 
lead to the bound Theorem E] (i). Most of the work in ^2] is in the proof 
of the bound (ii), which involves ruling out the possibility of other types 
of nucleation center with substantially higher probabilities than the event 
described above. 

The following set-up will be convenient. Let X be a random subset of 
in which each site is independently included with probability s. More 
formally, denote by Pg the product probability measure with parameter s 
on the product cr-algebra of {0, 1}^ , and define the random variable X by 
X{uj) = {x ^1? : uj{x) = 1} for oo G {0, A site x G ZMs said to be 
occupied if x G X. 

Consider the threshold growth model as in Theorem EJ For a set of sites 
K, let (K) denote the set of eventually active sites if we start with K as the 
set of initially active sites. We say that a set C is internally spanned 
if {X n K) = K. A rectangle is a set of sites of the form R{a,b]c,d) : = 
{a, . . . , c} X {6, . . . , d}; in particular we write R{c, d) = R{1, 1; c, d). Thus 

I{L,s) = Ps[R{L, L) is internally spanned). 
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Lower bound. Theorem E] (i) is a consequence of the following lemma, 
which is the analogue of Theorem 2 (i) in |12j . 

Lemma 11 

lim sup sup —s log J(m, s) < 2A. 

s— >0 m>l 

(Here A is as in Theorem Theorem 0] (i) is deduced from Lemma ITT] 
exactly as in the proof of Theorem 1 (i) in ^21, which in turn is based on an 
argument in [2]. The ideas follows; see ^2] or [2] for more details. Let 

S = S{L, s) be the event that R{L, L) contains an internally spanned square 
of side [s~^\ . The lemma implies easily that for L, s as in Theorem |3] (i) we 
have P5(>S') 1. On the other hand it may be shown that 

Ps(R{L, L) is internally spanned | S*) ^ 1, 

proving the required result. Next we turn to the proof of the lemma. 

Proof of Lemma [TH Let H = H{m) be the event that all the following 
occur: 

all sites in R{k, k) are occupied, 

the sites (m, 1), (l,m) are occupied, 

Ci, . . . , Cm-k-i has no /c-gaps, 

C'l,..., C'^-k-i has no /c-gaps, 

where Ci is the event that {k + i, j) is occupied for some 1 < j < i, and C[ 
is the event that (j, k + i) is occupied for some 1 < j < i. See Figure |2l 
If H occurs then i?(m, m) is internally spanned; to check this note that we 
may find an increasing sequence of internally spanned rectangles R{i,j) with 
i,j e [k,m] and \i — j\ < k, starting with R{k, k) and ending with R{m,m). 
On the other hand, we have 

Ps{H) > s^'+'PsiAkf 

where Ps{Ak) is as in Theorem |21 Therefore applying Theorem |21 yields the 
result. □ 
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Figure 2: An illustration of the event H with /c = 3, m = 12. 
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Upper bound. We now turn to the proof of Theorem E] (ii). We start 
by defining a new cellular automaton model called the enhanced model. We 
will explain how to prove the statement in Theorem |3] (ii) for the enhanced 
model, and this will imply the statement for the original model. 

For a finite set of sites F, let TZ{F) be the smallest rectangle containing 
F. Recall the definition of Nk in Theorem |3] The enhanced model is 
defined by the following rules. 

(i) A site that is active at time t remains active at time t + 1. 

(ii) For any site z, if F is any set of k sites in Nk{z) that are active at time 
t, then all sites in TZ{F) are active at time t + 1. 

(iii) For any site z', if F' is any set of two sites in N2{z') that are active at 
time t, then z' is active at time t + 1. 

(iv) Otherwise, inactive sites remain inactive. 

We define (■) and internally spanned for the enhanced model in the same 
way as for the original model. It is clear that any set that is internally 
spanned for the original model is internally spanned for the enhanced model. 
Hence it is sufficient to prove that the statement in Theorem |3] (ii) holds 
for the enhanced model. On the other hand, the enhanced model has the 
following useful property not shared by the original model: if ii' is a connected 
set of sites (in the sense of nearest neighbor connectivity on Z^) then (K) = 
TZ{K). Furthermore, the rectangle TZ{F) in (ii) is of course connected. 

We work with the enhanced model from now on, and our goal is to prove 
the statement in Theorem |^ (ii). (As a consequence, both parts of Theorem 
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Elhold for both the enhanced model and the original model, with the same A. 
This is in contrast with the situation for critical probabilities of percolation 
models, where any "essential enhancement" of a model strictly lowers the 
critical probability; see P). The proof follows the lines of ^2); and uses the 
following four lemmas, which we prove below. The first three are determin- 
istic properties of the enhanced model, and are analogues of Lemmas 9,10 
and Proposition 30 in ^2] respectively. 

Lemma 12 For the enhanced model, if a rectangle is internally spanned then 
it has no k consecutive unoccupied columns and no k consecutive unoccupied 
rows. 

The long side of a rectangle R = R{a, b; c, d) is defined as long(i?) = 
max{c — a + l,d — b+1}. 

Lemma 13 Let i be a positive integer and let R be a rectangle satisfying 
long(i?) > i. For the enhanced model, if R is internally spanned then there 
exists an internally spanned rectangle R' <^ R with long(i?') G [£,2i + 2k]. 

For rectangles R, R' we write dist(-R, R') = min{||x — x'\\oo '■ x & R, x' & 
R'}. 

Lemma 14 Let R be a rectangle with \R\ > 2. For the enhanced model, 
if R is internally spanned then there exist distinct non-empty rectangles 
Ri, . . . , Rr, such that 

(i) 2<r <k 

(a) the strict inclusions Ri <Z R hold for each i, 
(m) {RiU---URr) = n{Ri U---URr) = R, 

(iv) dist{Ri, Rj) < 2k for each i,j, 

(v) Ri, . . . ,Rr are disjointly internally spanned; that is, there exist dis- 
joint sets of occupied sites Ki, . . . , Kr with (Ki) = Ri for each i. 

Lemma 15 Fork > 1, the function fk,k+i is continuously differentiable and 
concave. 

Once Lemmas IT^HT^ are established, the proof of Theorem E] (ii) for the 
enhanced model requires only minor modifications to the proof in ^2]; so 



22 



we omit the details. The main modifications follows, referring to the 

terminology in that article. The definition of a rectangle being "horizontally 
(respectively vertically) traversable" should be replaced with statement that 
no k consecutive columns (respectively rows) are unoccupied. The function 
"(7(2;)" in [12] should be replaced with — log /fc,fc+i(e~^) in our notation (as in 
the proof of Theorem|2I). The variational principles in Section 6 of ^21 require 
that this function be convex; this follows from Lemma E| Finally, in the 
definition of a "hierarchy", a vertex may have up to k children (corresponding 
to having up to k rectangles in Lemma IT^ above) . and a vertex is declared a 
"splitter" if it has two or more children. 
Finally we give proofs of Lemmas IT^HT^ 

Proof of Lemma I12L Suppose that a rectangle R has k consecutive 
unoccupied columns, say. Then it is easy to see that even if all other sites in 
R is occupied, then no site in these columns lies in (R). □ 

Proof of Lemma I13L This is a straightforward corollary of Lemma IT^ 
Apply Lemma El to R, choose the resulting rectangle Ri with the longest 
long side, then apply the lemma again to Ri, and so on, until the required 
R' is obtained. □ 

Proof of Lemma [TH The proof is adapted from [12]. Let K = Rn X. 
By assumption we have (K) = R. We will construct (K) via the following 
algorithm. We introduce a new time parameter T not directly related to t 
in the definition of the model. For each time step T = 0, 1, . . . , r, we shall 
construct a collection of niT rectangles i?f , . . . , R'^j., and corresponding sets 
of sites Kf, . . . , -ft'j^j,, with the following properties: 

(i) Kf, . . . , are pairwise disjoint; 

(ii) Kr C K; 

R RT = {Kry, 

(iv) if z 7^ J then Rf ^ Rj; 
(v) K CU^ C (K), where 

rriT 
i=l 
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Initially, the rectangles and sets of sites are just the individual sites of K. 
That is, let K be enumerated as K = {xi, . . . and set mg = n and 

= Kf = {xi}, so that in particular 



The final set of rectangles will have the property that 

= (K) = R. (17) 

We call a set of distinct rectangles {-Ri, . . . , Rr} a clique if they satisfy 
conditions (i) and (iv) in Lemma El and in addition (Ri U ■ ■ ■ U Rr) = 
lZ{Ri U ■ ■ ■ U Rr). The algorithm proceeds as follows. Suppose i?f , . . . , R^^ 
and A'f , . . . , K"^^ have already been constructed. 

Step (I). If there does not exist a clique among i?f , . . . , -R^^, then stop, and 
set T = T. 

Step (II). Suppose there does exist a clique; after reordering indices if neces- 
sary, denote it i?f , . . . R^. Write R' = {Rju- ■ -URj) = 7^(i^lU■ ■ -URr), 
and K' = KfU---U Kj. 

Step (III). Construct the state (i2f+\ ),..., (i?^+^^, J at time 

T + 1 as follows. From the list (_Rf , K^), . . . , (-R^t' -^my) at time T, 
delete every pair (_Rf , Kf) for which R^ C R! . This includes (-Rf , Kj) 
for i = 1, . . . , r, and may include others. Then add (-R', K') to the list. 

Step (IV). Increase T by 1 and return to Step (I). 

It is straightforward to see that properties (i)-(v) are preserved by this 
procedure. Also mrp is strictly decreasing with T, so the algorithm must 
stop eventually. We must check that (fTTj) is satisfied. Suppose not; then 
there exists a site m R\U'^ which would become active in one step of the 
enhanced model if we start with active. Hence there exist z, F or 2', F' as 
in rules (ii),(iii) of the enhanced model. But the sites of F or F' must lie in 
at most k different rectangles . . . , i?^ say; then it is easy to see that these 
rectangles form a clique, so the algorithm should not have stopped. We claim 
that furthermore we must have ttIt- = 1 and R\ = R. If not, since LC^ = R, 
there must exist two distinct rectangles R^, i?J whose union is connected, but 
these two form a clique, so again the algorithm should not have stopped. 

Finally, considering the last time step of the algorithm (from time r — 1 
to time r) we obtain a set of rectangles with all the required properties. □ 



24 



Proof of Lemma fTHl It is sufficient to check the following 

(i) f"{x) < for X ^ p, 

(ii) /' is continuous at p. 

Let 4>{x) = x''{l — x), so that = and recall that (j)'{x) = only 

ai X = p = k/{k -\- 1). Differentiating (f){x) = (f)[f{x)] twice gives for x ^ p: 

fix) = '^^^^^ (18) 
^ r(x)-r[/(x)][r(a:)]^ 



where 



•0(0;) 



(^"(x) k {k-l)-{k + l)x 



[(f)'{x)f x^ [k-{k + l)x]2 ■ 

We have < for x < p and 0'(/(x)) > for x > p, therefore to 

prove (i) we need to show that 

il^ix) > tplfix)] for < X < p (19) 

and 

ip(x) < ip[f{x)] for p < X < 1. 

Since f[f{x)] = x, these two statements are equivalent, so we will prove the 
ffist. 

The ffist observation is that ip is decreasing on [0, p) and increasing on 
(p, 1] . To see this, compute 

k k-l + {k + l)[{k + l)x-{k-l)Y 



ip'{x) 



x'^+i [x(A; + 1) - 



and note that the numerator is nonnegative for all x and the denominator is 
negative for x < p and positive for x > p. This immediately gives (fT^ for 
0<x< (A;-1)/(A; + 1) < p: 

Hx) > ^— ^ ) =0> -2k = ^{1) > ij[f{x)]. 



k+1 
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To prove (|19|) for {k — l)/{k + 1) < x < p, we need the following two 
facts: 

<P{p - cr) > (p{p + cr) and ?/'(p - a) > ij{p + a) (20) 

for 0<cr<l/(A;+l). Before checking this, we will show that they imply 
(0 for {k - l)/{k + l) < X < p. Let X = p - a, so that < a < l/{k + 1). 
Then /(x) < p + cr = 2p — x by the first statement in ((201) • Using the second 
statement and the fact that ip is increasing on (p, 1] then gives 

^(x) = ^(p -a)> ^{p + a)> ^[/(x)]. 

It remains to check (|2Up. A little algebra shows that 

^P{p -a)- ^(p + a) = ^^Z^k ' 

so the two statements in ()20|) are equivalent. To check the first, compute 
d 



da 



p - a) - 0(p + a)] = a{k + 1) [(p + af-^ - (p - af~^] , 



so that 0(p — a) — 0(p + a) is increasing in a. Since this quantity is zero 
when 0" = 0, it follows that it is nonnegative for < a < l/{k + 1). 

We have established (i). This implies that the limits f'{p+), f'{p—) exist 
(but are possibly infinite); to check (ii) it remains to show that they are 
equal. Since (p" is continuous and non-zero at p, applying THopital's rule to 
fITHll gives 

f'ip-) - "^"^'-^ - ^— (21) 
^"r(p+)/'(p-)"/'(p-)' ^ ^ 

And similarly for f'{p+). But / is decreasing, so we must have f'{p—) = 
/'(p+) = -1. □ 
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